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Abstract 

The strongly coupled phase of Yang-Mills plasma with arbitrary gauge group is studied in a 
T-matrix approach. The existence of lowest-lying glueballs, interpreted as bound states of two 
transverse gluons (quasi-particles in a many-body set up), is analyzed in a non-perturbative scat- 
tering formalism with the input of lattice-QCD static potentials. Glueballs are actually found to 



■ be bound up to 1.3 Tc- Starting from the T-matrix, the plasma equation of state is computed 

by resorting to Dashen, Ma and Bernstein's formulation of statistical mechanics and favorably 
compared to quenched lattice data. Special emphasis is put on SU(iV) gauge groups, for which 
analytical results can be obtained in the large- limit, and predictions for a G2 gauge group are 
also given within this work. 
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I. INTRODUCTION 

More than two decades after pioneering works [l, 2], the phenomenology related to the 
deconfined phase of QCD, i.e. the quark-gluon plasma (QGP) is still a fascinating topic 
both experimentally and theoretically. On the experimental side, the QCD matter was or is 
studied in heavy-ion collisions (RHIC, SPS, FAIR, LHC). These experiments seem to show 
that the QGP behaves like a perfect fluid. On the theoretical side, the study of QCD at 
finite temperature deserves also much interest because it is a challenging problem in itself 
and because of the many connections with experiments. 

The aim of this work is to study the thermodynamic features of QGP by resorting to 
a T-matrix approach. The power of this approach is that the bound states and scattering 
states of the system can be studied in a whole picture. Such an approach has already 
proved to give relevant results in the study of hadronic matter above the critical temperature 
of deconfinement (Tc) jsj] but has not yet been applied to compute the equation of state 
(EoS). This observable will be performed here thanks to the Dashen, Ma and Bernstein's 



formulation of statistical mechanics in terms of the 5- matrix (or T-matrix) {4]. Such a 
formulation is particularly well suited for systems whose microscopic constituents behave 
according to relativistic quantum mechanics. The QGP is indeed identified to a quantum 
gas of gluons and quarks, which are seen as the effective degrees of freedom propagating in 
the plasma. This assumption is actually common to all the so-called quasiparticle approaches 
jsl , with the crucial difference that the use of a T-matrix formulation allows us to investigate 
the behavior of the QGP in a temperature range where it is strongly interacting. This strong 
interaction means here that bound states are expected to still survive above Tc. 

Although the above formulation can be applied to the full QGP, this paper is dedicated 
to the description of the gluon plasma. Dealing with only one particle species simplifies 
drastically the problem while the main feature of the description, i.e. the explicit inclusion 
of interactions in a quasiparticle approach, is kept. Moreover, the pure gauge thermodynamic 
features (in particular, the EoS) are well-known in lattice QCD; This will allow an accurate 
comparison between our phenomenological approach and the lattice QCD calculations. 

A particularity of this paper is the generalization of the formalism to any gauge groups, 
with a particular attention for SU(A^) and the large- limit, and for G2. This group has 
originally attracted attention because, the center of G2 being trivial, models relating decon- 



finement to the breaking of a center of symmetry are no longer valid as for SU(A^). However, 
it still exhibits a first-order phase transition as SU(A^) does j^. Hence, G2 appears quite 
attractive from a theoretical point of view. 

The paper is organized as follows. Sec. II is dedicated to the presentation of the general 
quasiparticle approach based on the T-matrix formalism proposed in J|. In Sec. HI, the 
model is particularized to a Yang-Mills plasma with the inclusion of 2-body interactions 
and, in Sec. IV, useful analytic comments concerning the thermodynamic observables in 
the SU(A^) and G2 cases are discussed. The model parameters are fixed in Sec. V and the 
existence of the bound states inside the gluon plasma is discussed in Sec. VI. In Sec. VII, 
the computation of the EoS is presented. Finally, Sec. VIII is devoted to the conclusions 
and perspectives. 

II. T-MATRIX FORMALISM 
A. Generalities 



a 



The results of [J] can be summarized as follows: The grand potential Q, expressed as an 
energy density, of an interacting particle gas is given by (in units where h = c = ks = 



(1) 



In the above equation, the first term, Qq, is the grand potential of the free relativistic 
particles, i.e. the remaining part of the grand potential if the interactions are turned off. 
The second term accounts for interactions in the plasma and is a sum running on all the 
species, the number of particles included, and the quantum numbers necessary to fix a 
channel. The set of all these channels is generically denoted u. The vectors fl = (/xi, fi2, • • • ) 
and = {Ni, N2, . . .) contain the chemical potentials and the particle number of each 
species taking part in a given scattering channel. 

The contributions above and below the thresholcj^ M^, are separated. Below the threshold, 
one has Q^, the grand potential coming from bound states, seen as free additional species 
in the plasma and appearing as poles of the S-matrix. Above the threshold, one has the 



^ Within this approach, the threshold is the summation on the mass of aU the particles included in a given 
channel v. 
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scattering contribution, where the trace is taken in the center of mass frame of the channel v 
and where S is the S-matrix, depending in particular on the total energy e. The symmetrizer 
S enforces the Pauli principle when a channel involving identical particles is considered, and 
the subscript c means that only the connected scattering diagrams are taken into account. 
Notice that K2{x) is the modified Bessel function of the second kind, that /3 is linked to the 
temperature T thanks to /3 = 1/T, and that the notation Ad^B = A{dxB) — {dxA)B is 
used. 

By definition, 5" = 1 — 27ri 6{e — Hq) T, where T is the off-shell T-matrix and where 
Hq is the free Hamiltonian of the system. A convenient way to compute T is to solve the 
Lippmann-Schwinger equation for the off-shell T-matrix, schematically given by 

r = v + vGor, (2) 

with Go the free propagator and V the interaction potential. 

Once the T-matrix is known, two problems can be simultaneously addressed: The exis- 
tence of bound states in the plasma and its EoS. The T-matrix formalism has the advantage 
of treating bound and scattering states on the same footing, and is particularly suited for 
the present situation where we expect bound states to become less and less bound when the 
temperature increases, eventually crossing over and melting into the continuum. This dis- 
sociation mechanism has been shown to provide considerable threshold enhancement effects 
in heavy quark anti-quark correlation functions jsl. 

The plasma EoS is obtained from ([1]). Then, the pressure is simply given by 

P = (3) 

and the other thermodynamic observables can derived from p. For example, the trace 
anomaly ( A = e — 3 p) and the entropy density (s) read 

where V is the volume of the system. For later convenience, the thermodynamic quantities 
will be normalized to the Stefan-Boltzmann pressure, which is defined as 

PsB = - lim ^lo, (5) 
rrii being the masses of the particles propagating in the medium. 
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B. Interaction potential 



The explicit computation of fl obviously requires the knowledge of the on-shell T-matrix 
that can be derived in particular from A key ingredient of the present approach is 
thus the potential V, encoding the interactions between the particles in the plasma. In 
the following, V is chosen as pairwise: For a n-body channel, V = Yli<j ^ij with Vij the 
potential between two particle species i and j. Having in mind the building of an effective 
framework describing the deconfined phase of a non-abelian gauge theory, each particle 
composing the plasma should be in a given representation of the considered gauge (or color) 
group. It is therefore reasonable to assume that the potential V between two particles in 
the representations Ri and Rj of the considered gauge group has the color- dependence of a 
(screened) one-gluon-exchange potential, that is, in momentum space, 

V^, = MR^-MR^asviP,q,q), (6) 

where Mr denotes the generator of the considered gauge algebra in the representation R, 
and where the real function v only depends on the temperature and two momenta (no 
dependence on the mass or other attributes of the particle). We keep the name gluon for the 
gauge particle even if the gauge group can formally be arbitrary. In the above definition, 
it has to be remembered that as = (7^/47? and that = A/Cg''"', adj being the adjoint 
representation of the gauge group under study and C^f being the value of the quadratic 
Casimir in the representation R. Note that in the case of SU(A^), A is the 't Hooft coupling 
(fixed in the large- limit). 

Introducing quadratic Casimirs, one can rewrite ([6]) as 

Vij = asv^ Kc-ij V, (7) 

with C the pair representation and 

- ^-^j • (») 

Again, the real function v = v{(3, g, q') only depends on the temperature and on two mo- 
menta - an explicit form for v will be given later. The validity of the form (I7j) for Vij has 
partially been checked in pure gauge SU(3) lattice calculations, showing that the static po- 
tential between two sources, in different representations and bound in a color singlet, follows 



the Casimir scaling expected from a process of one-gluon-exchange type [7|]. The pecuhar 
scahng ([7]) also leads to a relevant large- iV behavior of the EoS when the gauge group SU(A^) 
is chosen, as it will be shown in Sec. IIVI 

Among the various possible representations, the case where a singlet (denoted •) appears 
in the tensor product Ri ^ Rj is particularly relevant: Since C* = and the other quadratic 
Casimirs are positive, the singlet is the most attractive channel in any two-body scattering 
process, so the most favorable one for the formation of bound states. Such two-particle 
bound states should presumably be the lowest-lying ones and, being color-singlet, would 
give rise to low-lying glueballs or mesons for instance. 



C. Born approximation 

The scattering term in ([T]), given by 

(9) 

can be considerably simplified by using the Born approximation, i.e. by noticing that if 
the interactions are weak enough, T = V + O(y^). Such conditions are generally expected 
to be valid at high enough temperatures, where the typical interaction energy is small 
with respect to the typical thermal energy of the particles. Note also that, in some cases, 
this approximation can be relevant when the factor kcaj is negligible, irrespective of the 
temperature: Such cases will be encountered when the gauge group is SU(A^) (see Sec. HV|) . 
To the first order in V, (jH]) becomes 

^s = Y.Vv^2 dee'K,i(3e)TT^dM^-Ho)V)l + 0{V'). (10) 

Let us write explicitly u = (n,!)), where n is the total number of particles involved in 
a given scattering channel, and where u are the remaining quantum numbers. A useful 
remark to be done at this stage is that the pairwise structure of V causes V\^ to be always 
vanishing excepted in two-body channels. Here, at the Born approximation, n is always 
equal to 2. Then, Tr;^ de{S{e — Ho)V)\^ = Tr/> nc,ij de{Sv), with 6 = 6{e — eij{q)) and €ij{q) = 
\/ + + \J 'f' + ^"j- Note that the color channel C and the particles species i, j are part 
of u. After having extracted from the trace the color and J^*" dependences ( if the charge 
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conjugation is not relevant), one is led to 
^^ = E^-^E(2^+l)E^™'^"^c,.W dee' K2mTT,d,{6vjPc) + 0{V^), (11) 

with dim C is the pair representation dimension, Tr^ the remaining trace on the momentum 
space and Vjpc the potential with the angular symmetry of the considered channel. 

Among the various summations to be performed in two are of particular interest: 
The one over the different interacting species, that can be denoted ^"^^ 
over the representations appearing in Ri ® Rj, that is Because of i^c,ij, dUD is thus 

proportional to a factor ^^^dimC Kc,ij ^oi a given pair i, j in a given J^"^ channel. When 
the combinations of species does not have to respect a symmetry principle, this last sum 
runs on all the representations appearing in Ri Rj', one can then show that 

^dimC Kc,ij = 0. (12) 
c 

Indeed, it is known in group theory that the second order Dynkin indices in a tensor prod- 
uct obey a sum rule that can be rewritten using our notations as I^' dim Rj + 1^^ dim Ri = 
J'' 81]. Using C2 = {dimadj /dimR)I^ 8|], one straightforwardly shows that ( fT2|) holds. 
Note that (fT2l) and ( fTTl) are thus a priori nonzero when a symmetry principle has to be 
respected: The summation cannot then be performed on all possible color representations. 



III. YANG-MILLS PLASMA 
A. Grand canonical potential 

Let us now particularize the general formalism presented in the previous section to a 
genuine Yang-Mills plasma, i.e. with no matter fields. The bosonic degrees of freedom 
propagating in the plasma are then quasigluons, that are transverse spin-1 bosons in the 
adjoint representation of the gauge group. The baryonic potential can be set equal to and 
according to standard formulas in statistical mechanics, one has 

Qo = 2 dim adj uo{mg), (13) 

where the quasigluons are a priori supposed to have a mass nig, and where 

cuo(m) = ^J^"^dk P In (1 - e-^v^^j (14) 
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is the grand potential per degree of freedom associated to a bosonic species with mass m. 
Equation ([5]) leads to 



PsB = ^^dim adj. (15) 

Let us recall that in the following, the term gluon indifferently denotes the gauge field of 
Yang-Mills theory and the quasigluons. 

The sum appearing in ([1]) now explicitly reads J2jpc: where Ug is the number 

of gluons involved in the interaction process. As soon as > 2, the determination of the 
allowed color channels and of the correct symmetrized gluon states generally becomes a 
painful task, to which the problem of finding the T-matrix in many-body scattering must 
be added. Intuitively, one can nevertheless expect the dominant scattering processes to be 
two-gluon ones, and thus only consider = 2 in a first approach. After simplification, the 
grand potential ([1]) eventually reads 



= 2 dim adj ^0(^3) + dimC (2 J + 1) I 

c jp I 



uj{Mc,jp) (16) 



27r2/3' 



1 f°° r 

— / dee^ K2{f3e) Tic,jp {5ReT)' - 27r((5Rer)(5Imr)' - ((5Imr)(5Rer)') 



where the symbol "prime" is the derivative respective to the energy and Mq jp is the mass of 
the two-gluon bound state with color C and quantum numbers , if it exists. The index C 
in the J^'" channel is dropped since the charge conjugation is always positive for a two-gluon 



state jo]. In the remaining trace, it is understood that the T-matrix has been computed in a 
given two-body channel with color C and quantum numbers J^, and that the Dirac 5 reads 
5{e — 2e(g)), with the dispersion relation e(g) = ^Jq'^ + m'^. Note also that, in connection 
with nuclear many-body approaches, ( fT6|) can be rewritten in terms of a weighted thermal 
average of scattering phase shifts by means of unitarity of the S'-matrix. The computation 
of the two-gluon T-matrix is explained in detail in the following section. 



B. Helicity states and the Lippman-Schwinger equation 

1. Two gluon states 



Jacob and Wick's helicity formalism 10| can be applied to describe a two-gluon state, 
where the gluons are seen as transverse spin-1 particles. Let us generically define \ip{p, ^)) = 

8 



a\{p) |0) the quantum state of a particle with momentum p, spin s, and hehcity A. If the 
particle is transverse, only A = ±s is allowed, while all the projections from —s to +s are 
allowed if the particle has a usual spin degree of freedom. Then it can be deduced from [K 
that the quantum state 



|j^,M;Ai,A2,r/) 



1 

TIL 



|J,M;Ai,A2)+r/|J,M;-Ai,-A2) 



(17) 



with rj = ±1 and 

|J,M;Ai,A2) 



'2J+ 1" 


2 /■27r 




Ait 


/ 

Jo 


^0 



de sine vii,^_,^{<f>,i 

xal{p)ali-p) |0) 



(18) 



is a two-particle helicity state in the rest frame of the system which is also an eigenstate of the 
total spin J and of the parity, i.e. = J{J + 1), Jz = M, and P = e?7i?72(— 1)"'~*^~^^ with 
rji and Sj the intrinsic parity and spin of particle i. Moreover, J > |Ai — A2I. In the above 
definition, R{a,f3,'y) is the rotation operator of Euler angles {a,/3,7} and V'lj ^^{a, f3,'y) 
are the Wigner Z)-matrices. The coordinates {6, 0} are the polar angles of p. When both 
particles have a spin degree of freedom, the helicity basis, spanned by the helicity states (ITTI) . 
is equivalent to a standard |^'^"'"^Lj) basis up to an orthogonal transformation Hj]. When at 
least one of the particles is transverse, both basis are no longer equivalent, but the helicity 
states can still be expressed as particular linear combinations of states 10|. This 

will be convenient in view of future computations. 



When the two particles are identical (mi = m2 = m, s 



the action of the permutation operator P12. One finds 10 1 



S2 = s), it is relevant to study 



[1 + (-1)2-Pi2] I JP, M- Ai, A2, Ti) = I J^, M; Ai, A2, ri) + (-1)^ | J^, M; A2, Ai, r/) , (19) 

where the operator [1 + (— l)^*Pi2] = iS is nothing else than a projector on the symmetric (s 
integer) or antisymmetric (s half-integer) part of the helicity state. It is readily seen in (fT9|) 
that symmetrizing the state will eventually lead to selection rules for J (this is particularly 
clear if one sets Ai = A2). When extra degrees of freedom are added, it is also of interest to 
use the antisymmetrizer [1 — (— l)^^Pi2] = .4, as done in Table [B 



A general discussion about the two-gluon helicity states can be found in 12|, to which 
we refer the interested reader. For the present work, it is sufficient to recall that four 
families of helicity states can be found, separated in helicity singlets \S±] J^) and doublets 



TABLE I. Symmetrized and antisymmetrized two-gluon helicity states, following the notation 
of [l^. with the corresponding quantum numbers and averaged squared orbital angular mo- 
mentum. 



State 


Symmetrized Antisymmetrized 






\s+ 


JO 


(even-J > 0)+ (odd-J > 1)^ 


J(J + 1) 


+ 2 






(even-J > 0)" (odd-J > 1)+ 


J(J + 1) 


+ 2 






(even-J > 2)+ (odd-J > 3)" 


J(J + 1) 


- 2 






(odd-J > 3)+ (even-J > 2)" 


J(J + 1) 


- 2 



\D±; J'^) following the pioneering work 13 1. The corresponding quantum numbers are given 
in Table [H as well as the average value of the squared orbital angular momentum, {p^'^i 
computed with these states. 

The averaged orbital angular momentum is an interesting quantity since it helps to glob- 



ally understand the mass hierarchy of the glueball spectrum [12|. Moreover, in a naive 
nonrelativistic picture, it estimates the strength of the orbital barrier in scattering theory. 
For obvious numerical reasons, all the possible channels contributing to can not be 
included, that is why it is of interest to find the channels that will presumably contribute 
the most, i.e. those with the lowest value of [L'^)- First, one has the symmetric states 



\s-o~) 



1/2 



1/2 



'Do) 



(20) 
(21) 



expressed in a standard \'^^~^^Lj) basis, with (^-^^^ = 2. In the singlet channel, they corre- 
spond to the 0+"'' and 0""'' glueballs respectively, which are indeed found to be among the 
lightest ones at zero temperature, see e.g. the review [l^. Then, with ( L'^\ = 4, one has 
the symmetric state 



1/2 



1/2 



' 1 ■ 

35 



1/2 



(22) 
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and the antisymmetric states 





"2" 

3 


1/2 






'1" 

3 


1/2 





15 

2 

3 



1/2 



1/2 



\D_;2-) 



1/2 



^^2> 



rini/2 



'Pl>^ 
^/^i>, 



1/2 



^^1> 



(23) 
(24) 
(25) 



The above three states do not exist in the singlet channel, but the symmetric state 
corresponds to the 2+"^ glueball in the singlet channel; the O^"*" and 2~^~^ are indeed the 

n 

lightest states at zero temperature p^. Only the color-symmetric channels with the lowest 
value of ^-^^^ will be kept in the following study, which aims at being a first step toward a 
description of the Yang- Mills plasma within a T-matrix formulation. 



2. Lippman-Schwinger equation 

Solving ([2]) is a crucial technical part of this work since it eventually leads to the on-shell 
T-matrix. As it will be discussed in SecJVl the potential to be used is known in position 
space and has firstly to be Fourier-transformed. For a potential with spherical symmetry in 
configuration space, we use 

r°° sinfOr) 

V{q,q,eq,q>)=A-K I drrV{r) — where Q = ^ q^ + q'^ - 2qq' cos 6 q,q' , (26) 

Jo Q 

and where Oq^qi is the angle between the momenta q and q' . 

Since two-gluon interactions are considered, the basis states are two-gluon helicity states, 
given in the above section. As we assume V to be spin independent (see Appendix [B]), only 
the orbital angular momentum containing of the helicity states has to be taken into account. 
According to a standard integration, the L-wave part of potential (126|) reads 

Fi(g,g') = 27r J ^dxPLix)V{q,q',x), (27) 

where Pl is the Legendre polynomial of order L and x = cos 6q q/. Our choice is to focus 
on the scalar, pseudoscalar and tensor scattering channels, for which one can compute from 
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([20])-(l22]) that 

= lVo{q,q') + ^V2{q,q'), (28) 
Vo~{Q,q') = Vi{q,q'), (29) 

^2+(g,g') = lvo{q,q') + ^\/2(g,g ) + ^^4(9, g')- (30) 

Note that once Vjp{q,q') is known, that is the potential in a given scattering channel, 
the off-shell T- matrix can be computed from ([2]) as follows in 3|: 

T{E;q,q') = Vjp{q,q') + — / dk ^ Vjp{q, k) Go{E; k) T{E; k, q'), (31) 

where the two-gluon propagator reads 

1 

= ^EV4-e(A;)2-2.e(A;)E, ^^^^ 



with the gluon dispersion relation e{k) = ^^Jk'^ + rn?g. Note that the normalization conven- 
tions of the T- matrix are not the same as the ones in (see Appendix O). The parameter 
S/ accounts for the imaginary part of the gluon self-interaction, whereas the real part is 
reabsorbed in the effective gluon mass. A more complete calculation of the gluon self-energy 
would require summing the T-matrix over the gluon thermal distribution self-consistently 
in a Brueckner-Hartree-Fock scheme (schematically, Ilg = j f^TDg with the gluon dis- 
tribution and Dg the gluon single-particle propagator). We leave such determination for a 
future work. In the present evaluation, we shall approximate the gluon self-interaction by 
using an effective, in-medium gluon mass (to be discussed in Sec. |V]) together with a small 
imaginary part for numerical purposes (we use S/ = 0.01 GeV as in j^). 

Once T{E; q, q') is known, the on-shell T-matrix is readily obtained as T{E; qE, qs), with 



qE = \J — m?g. The Haftel-Tabakin algorithm is used to solve ( 13T|) |l5| . The momentum 
integral is discretized within an appropriate quadrature, thus turning the integral equation 
into a matrix equation, namely, 'Y^J^ikTkj = Vij, where, schematically, J-" = 1 — wVG (and 
w denotes the integration weight). The solution follows trivially by matrix inversion. It can 
be shown that the determinant of the transition function J-" (referred to as the Fredholm 
determinant) vanishes at the bound state energies, which provides a numerical criterion 
for solving the bound state problem. This strategy has already been successfully used to 
compute T- matrices in the case of quark-antiquark scattering j3|. 
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IV. THERMODYNAMIC OBSERVABLES WITH SU(iV) AND G2 



A. Pure gauge sector 

1. SU(N) case 

The explicit computation of VL^'^\ given by f|T6l) . obviously requires the knowledge of 
the on-shell T-matrix, that can be derived in particular from ([2]). In this last equation, 
Go is the propagator of two gluons, that has been discussed in Sec. IIIII It is only worth 
saying that Go = 0(1) with respect to the number of colors since nig is assumed to be 
0(1) (see Sec. IV AD . The color dependence of the T-matrix actually comes from the two- 
gluon interaction potential only. More precisely, the color- dependence of the potential is all 
included in the factor ([8]), reading in the present case 

I^C;gg = ^ • (33) 

The subscript gg is used to recall that two-gluon interactions are concerned in the above 
formula, and G| = ^2'^"' = in the SU(A^) case. 

The adjoint representation of SU(A^), to which gluons belong, can be written as the 
(A^ — l)-component vector (1,0,..., 0, 1) in a highest weight representation, corresponding 
to a Young diagram with 1 column of length — 1 and 1 column of length 1. More 
generally, (ai, . . . , a^, . . . , aN~i) corresponds to a Young diagram with columns of length 
k. The tensor product of the adjoint representation by itself gives the allowed two-gluon 
color channels: 

(1,0,. ..,0,1)0(1,0,. ..,0,1) = 

•^©(1,0,..., 0, 1)^ © (2, 0, . . . , 0, 2)^ iV > 2 

© (1, 0, . . . , 0, 1)^ © (0, 1, 0, . . . , 0, 2)^ © (2, 0, . . . , 1, 0)^ iV > 3 

©(0,1,0,..., 0,1,0)^ N>A. (34) 

The superscript S/A denotes a symmetric/antisymmetric channel. The first /second/ third 
line exists as soon as N > 2/3/4. Note that in the special case N = 2, the above tensor 
product reduces to (2) © (2) = (0)'^ © (2)^ © (4)'^, and one recovers usual spin-coupling rules. 
The dimensions and color factors of the representations appearing in ( 134|) can be found in 
Table Ml 
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In the singlet channel, one has k, = — 1 for any A^. It is such that T = 0(1) since 
V = 0(1). Consequently, the properties of glueballs in singlet above the deconfinement 



temperature are not dependent of TV, in agreement with [16| , where it is suggested that this 
argument is even gauge-group independent. The singlet finally brings a contribution 0(1) 
to fi'-^^ since its dimension is 1. 

Using the same arguments as for the singlet, one finds that the adjoint channels also lead 
to a T- matrix that is A^-independent. They may lead to bound states since the potential 
is attractive, though less strongly than for the singlet. The symmetric adjoint channel 
will presumably be the most favorable for the formation of bound states since it demands 
a completely symmetric spin-space wave function for the two-gluon state in virtue of the 
Pauli principle, and the most attractive channels are indeed symmetric. Note that this 
symmetric color channel is actually absent for N = 2. In the adjoint channel, T = 0(1) since 

V = 0(1) but, unlike the singlet, its contribution to fi*^^^ is 0(A^^) since dim(l, 0, . . . , 0, 1) = 
N'-l. 

The two remaining channels with nonzero potential, namely (2,0,..., 0, 2) (the 27 for 
SU(3)) and (0, 1, 0, . . . , 0, 1, 0) (only when > 3), have in common that they are symmetric 
and that their color factor scales in 1/A^, thus vanishes in the large- limit. The fact that 

V = 0{1/N) in both cases leads to the exact large- A^ result 

r =V + VGoV + 0{N-'^), (35) 

or 

r=±^v + ^^vGoV + OiN~'), (36) 

the lb coming from one channel or another. Because of the weakness of V at large- A^, one 
can reasonably suppose that even the attractive channel (0, 1, 0, . . . , 0, 1, 0) does not lead to 
the formation of bound states. For the two channels under consideration, 



o^^eRe(±^ + ^UO(Ar-). (37) 



"A^ A^2 

One sees in (1161) that the contributions of both channels have to be summed and, since they 
are symmetric, the sums on the allowed is identical in both cases. This causes the term 
in 1/A^ to vanish in the trace at large- A^ limit, the first nontrivial one being in l/N"^, leading 
to an overall contribution to fi'-^^ scaling as A^^ because the dimension of both channels scale 

as A^^ 
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Although the color singlet is relevant in view of studying glueballs, it does not bring any 
contribution to the EoS at large- A^. So, the large- iV EoS is dominated by free gluons and 
scattering processes above threshold in colored channels. The more is large, the more 
important is the gap between the confined phase and the deconfined one, whose EoS scales 
as A^^. It is indeed known that the large A^-case corresponds to a strongly first-order phase 
transition (A^ = 3 is already weakly first order) 17 1. 



2. G2 case 

Another interesting group under consideration is G2 which is also the best studied gauge 
group so far beyond SU(A^). The main features of this group are summarized in what follows. 

The adjoint representation of G2 has dimension 14, and reads (0, 1) in a highest weight 
representation. The two-gluon channels are then given by 

(0, 1) ® (0, 1) = + (0, 1)^ + (0, 2)^ + (2, 0)^ + (3, 0)^ (38) 

or, in terms of the dimensions, 14® 14 = 1 + 14 + 77' + 27 + 77. Using the same normalisation 



than in the SU(A^) case, the color factors respectively read [l8|, ll9[ Kc-gg = —1, —1/2, 1/4, 
—5/12, and 0. The color factors in the singlet and adjoint channels are equal to those of 
SU(A^), so the glueball properties are unchanged in the singlet and antisymmetric adjoint 
channels. The symmetric (2,0)"^ channel is almost as attractive as the adjoint one: It may 
lead to bound states. 



3. Scaling relations for SU(N) and G2 

Some interesting relations about the scaling of the EoS can be deduced thanks to the 
T-matrix. Let us write the on-shell T-matrix as T = Xlfc '^c gg where all do not depend 
on the color but rather on the other quantum numbers involved. The color dependence of 
the thermodynamic observables is then given by the quantities Ylc-A/s^^^^ '^C gg- Using 
the results of Appendix R] and Sec. IIV A 2t one can check that, for SU(A^) and G2, 
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^dimC, 



1 



3 



C:S 



dim adj, 



^ dim CggKl.gg = ^ dim adj, 



C:S 



4 



^dimCggK, 



3 

C;gg 



1 



dimarfj, 



C;S 



J2^^^<^gg4;gg 



C:A 



dimarfj. 



For SU(A^) at large- A^, the previous relations can be written 

k 



J2dunCgg4;gg = 



C:S 



j:dunCgg4;gg = N' (-0' 



Ok,l + ■^0k,2 



0(1). 



0(1] 



(39) 
(40) 
(41) 

(42) 



(43) 
(44) 



Again that means that the expected scaling like A^^ (actually like dim adj) of the EoS is 
found using the present approach. This can be viewed as a confirmation of the relevance of 
the chosen color scaling ([6]). 



B. Quarks and antiquarks in the 't Hooft limit 

Even if the rest of this study will be concerned with a genuine Yang-Mills plasma, it is 
worth making some comments about the possible inclusion of matter (quarks and antiquarks) 
in the model. Informations about the color channels appearing in interactions involving at 
least one (anti)quark are given in Tables IVl and IVTI remark that quarks (antiquarks) have 
been considered to be in the fundamental (conjugate) representation of SU(A^), as it is the 



case in 't Hooft large- A^ limit 20|. Other interesting large- A^ limits have been proposed, in 



which quarks belong to the two-index antisymmetric representation of SU(A^) for example 



2l|, but they will not be studied here. 

First of all, the quark-quark and antiquark-antiquark color factors are of order 1/A^: 
The Born approximation for the T-matrix becomes exact at large- A^ (at any temperature) 
and, since the dimension of any of the corresponding representations scales as A^^, the 
quark-quark and antiquark-antiquark interactions bring a term scaling as A^ to the grand 
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potential. More precisely, this term scales as Nj N at the Born approximation once the trace 
over the different ffavors is performed. It is shown in Sec. Ill CI that only the interaction of 
two identical species can contribute to Q in this limit. The number of quark flavors remains 
finite in the 't Hooft limit, which is the case under study here. 

The quark-antiquark interactions lead to a T-matrix which is 0(1) when the pair is in the 
singlet, or 0(1/A^^) when the pair is in the adjoint representation. In both cases however, 
the contribution to the grand potential scales as Nf{Nf + l)/2 at large- A^. So the quark- 
antiquark contributions to the thermodynamic observables is negligible with respect to the 
quark-quark and antiquark-antiquark ones in the 't Hooft limit. 

Finally, using similar arguments, one can show that the contribution of the quark-gluon 
and antiquark-gluon interactions to the grand potential scale as Nf N at large- A^. One 
concludes that, in the 't Hooft large- limit, the grand potential is dominated by the 
gluonic contributions only, scaling as A^^. 



V. PARAMETERS OF THE MODEL 



A. Potential and gluon mass 



Two ingredients are now missing to start numerical computations: The interaction po- 
tential between two gluons and the gluon mass. The procedure followed to fix the potential 
is similar to the one followed in 31] in the case of heavy quark-antiquark bound states. The 
first step is to take some input from lattice QCD, from which accurate computations of the 
static free energy of a quark-antiquark pair bound in a color singlet, Fi{r,T), are available. 



In particular, computations in quenched SU(3) lattice QCD can be found in [22]; they are 
especially relevant for our purpose since we focus on the pure Yang-Mills plasma. There is 
still debate on the proper potential term to use in phenomenological approaches, namely Fi 
or the internal energy Ui = Fi — TdrFi. An entropic contribution is subtracted from the 
free energy in Ui, causing the internal energy to be more attractive than the free energy, 
eventually leading to larger dissociation temperatures for bound states in the deconfined 
medium. Spectral function analysis of heavy quarkonia from lattice QCD simulations of 
Euclidean correlation functions typically suggest that the rjc and J/ip states may survive 
up to about 2Tc. Such values of the dissociation temperature can be accommodated if the 
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singlet internal energy is used in potential model calculations. [23|-|26|. As in [3|, the internal 
energy is used as potential term. The explicit expression of the internal energy Ui used in 
this work can be found in Appendix [Bl 

The assumed color scaling ([7]) allows to derive the two-gluon potential from the lattice 
quark- ant iquark one. Indeed, given Ui{r,T) computed in quenched SU^Niat) lattice data, 
the color factor of the singlet quark- ant iquark pair reads 

Ar2 - 1 

= -^k^- (45) 

According to (j7]), the potential (in position space) between two quasigluons in the color 
channel C is then given by 

V(r,T) = ^ [f/i(r,T) - U^{oo,T)] , (46) 

where the long-distance limit of the potential has be normalized to zero in order to ensure 
the convergence of the scattering equation and to perform the Fourier transform. This is 
actually a standard procedure in finite-temperature calculations. 



According to the suggestion made in [27|], the nonzero value of t/i(oo,T) should eventu- 
ally be responsible of an effective in-medium contribution to the gluon mass. The intuitive 
argument is that, when both gluons are infinitely separated, they no longer interact. There- 
fore, the remaining potential energy should be seen as a manifestation of self-energy effects 
induced by the surrounding medium. These effects are encoded in the model as a mass shift 
to the "bare" quasigluon mass, whose value has still to be fixed. Since t/i(oo,T) = 2mq(T), 
the adaptation to the gluon must be done by extracting the correct color-dependence. From 
HTL computations [28|], the self-energy color dependence is given by Cl/Cf at the first 
order when it is added in the propagator as a mass term (m^), that means here that 



- - ^^A(T). (47) 

So, mg(T) = 2A(T)/3 in the SU(3) case. 

In the same way as done for the two-body color scaling, A(T) is considered as universal 
and the gluon thermal mass reads 



W = W^A(T) = A(T), (48) 
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since d = Cf. So, S{T) is gauge-group independent. The effective in-medium gluon mass 
is finally given in our approach as 



5{Tf. 



(49) 



where the value mo has still to be fitted (see the following section). All the contributions 
are quadratically added as it is the case when one is dealing with bosonic propagators. The 
gluon mass is thus gauge-group independent. 

It is obvious that the problem of the gluon mass is far more complicated than the simple 
prescription f H^ . that has to be seen as valid in a first approximation only. A more re- 
fined gluon mass should probably be momentum-dependent. There is indeed an increasing 
amount of evidences favoring the existence of a dynamically generated gluon mass due to 
nonperturbative effects, at least at zero temperature. Such a dynamically generated gluon 
mass rriglp), with mJoo) = and ?Tig(0) finite, is favored by some lattice results in Lan- 



dau gauge, see e.g. 



29 



30|. Also nonperturbative field-theoretical calculations, using for 



example t 
theory 131 



le pinch technique, find a nonzero dynamically generated gluon mass in 3 + 1 YM 



32] . It is also worth quoting the recent Coulomb gauge study 



which is a first 



step in view of understanding the behavior of mg{p,T) at a nonperturbative level. From 
a different perspective, non-perturbative contributions to the gluon potential and mass are 



35| . Such 



analyzed at finite temperature in connection with the gluon condensates in 34 
improvements of the gluon mass are left for future works. 

The above discussion gives a more precise meaning to the term "quasigluon" used in this 
paper: It denotes transverse particles in the adjoint representation of SU(A^) that gain an 
effective mass mg{T) given by (l49l) and interact through the potential (146|1 . 



B. Zero temperature results with SU(3) 

Before performing finite-temperature computations, it is worth checking whether the 
values retained for the various parameters of our model may give relevant results at zero 
temperature or not. In particular, is the present T-matrix formalism able to reproduce at 



least qualitatively the features of the low- 
SU(3) lattice QCD at zero temperature ? 



ying glueball spectrum computed in pure gauge 



36| 



It is known that, at zero temperature and in quenched SU(3) lattice QCD, the potential 
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between a static quark-antiquark pair is compatible with the funnel form 37 1 



Vjir) = ar-^^. (50) 

In order to stay coherent with the potential above Tc, a = 0.141 (see Appendix [B]) and 
cr = 0.176 GeV^ (a standard value for the string tension) are used. The Fourier transform 
of Vf{r) is not defined: This flaw can be cured by making it saturate at some value Vgb, 
interpreted as a string-breaking value, that is the energy above which a light quark-antiquark 
pair can be created from the vacuum and break the QCD string. This scale is then subtracted 
and the potential effectively taken into account is V/(r) — Vsb, while Vsb is interpreted as an 
effective quark mass using the same arguments as those detailed in Sec. IV A[ According to 
the color scaling ([8]), the potential that should be used to describe the interactions between 
two gluons at zero temperature is 

Voir) = IVfir) - ^1 (51) 

when the gauge group is SU(3). In this case, the string breaking scale should rather be 
interpreted as the energy scale necessary to form two gluelumps, a gluelump being a gluon 
bound in the color field of a static adjoint source. It is known indeed that adjoint strin, 
breaking may be observed, and occurs at twice the lightest gluelump mass (~ 2 GeV) 



So, V^, = 2 GeV is used here, a value in agreement with lattice data showing that the mass 



of the lightest gluelump is given by 0.85(17) GeV [39 1. 

The only free parameter left to compute the T-matrix is the bare gluon mass, mo- Keeping 
the same structure as in Eq. (HUj) we have 



A3 \ 2 



m,(0)^ = mg+(%) , (52) 



2 

where again we have traded the subtracted potential at infinite separation distance (string 
breaking energy in this case) into a self-energy-like contribution to the quasi-particle gluon 
mass. We fix niQ = 0.7 GeV, which is a typical value for the zero-momentum limit of the 
gluon propagator at zero temperature. Advancing results, such a value will ensure both a 
correct agreement with the zero temperature lattice glueball spectrum (see Table HB . and 
an excellent agreement with the pure gauge EoS computed on the lattice (see Sec. IVIip . 

The results are given in Table HT] for channels discussed in Sec lIIIBl At least our model 
is able to reproduce the mass hierarchy of the lightest glueballs as well as the typical mass 
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TABLE II. Masses (in GeV) of the lowest-lying glueball states at zero temperature with the gauge 



group SU(3). Our results (third column), are compared to the lattice data of 



and to the Coulomb gauge QCD (CGQCD) study 



36( 1 (second column) 



40j] (last column). 



State 


Lattice [36] 


T-matrix CGQCD [40] 


0++ 


L73 (5)(8) 


2.17 1.98 


0-+ 


2.59 (4)(13) 


2.39 2.22 


2++ 


2.40 (2.5)(12) 


2.34 2.42 



scale of 2 GeV for those states. The accuracy of the model can be compared to Coulomb 
gauge QCD sharing formally many similarities with our T-matrix approach: The results 
of this last reference are also given in Table [ITl The agreement between lattice QCD and 
CGQCD is better, but it is worth saying that the parameters used in CGQCD have been 
chosen to reach an optimal agreement with the zero temperature lattice data, while here 
the values are mostly designed to give good results above Tc. So our model results are 
satisfactory in that sense. What makes us to find such a high mass for the scalar glueball 
is the quite small value a = 0.141 that has been taken (in order to fit static potentials 
above T^), while values as high as a = 0.4 have sometimes to be used to reach a good 
agreement between lattice data and effective approaches, see e.g. 4l|. The scalar glueball 
being dominantly a S-wave state, it is particularly sensitive to the strength of the Coulomb 
term and to the running of a with the temperature, that is neglected here. 

Finally, the extension of the above calculations to any gauge group is straightforward 
in our approach: The interested reader will find a discussion of such a generalization in 
19[, where it is shown that the lowest-lying glueball masses is gauge-group independent 
within a constituent framework. In particular, the lowest-lying glueball masses are found 



independent of in 19|, in agreement with what is observed on the lattice 42|. That is 
why the T-matrix masses given in Table |TT] are considered as valid for any gauge group too. 



C. Relevance of the parameters for other gauge groups 

At this stage, it is important to summarize the various parameters introduced and their 
possible dependence - or not - on the gauge group. Comparison with existing results when 
it is possible can also shed light on that issue. 
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First of all, the basic ingredient underlying the present study is the static quark-antiquark 
potential computed in finite-temperature quenched lattice QCD. The assumed one-gluon- 
exchange nature of the two-particle interactions leads to the universality of the momentum- 
dependent part of the potential, and to a well-defined prescription for its gauge-group de- 
pendence. Similarly, the gluon thermal mass has a peculiar color scaling originating in its 
interpretation self-energy term. 

More freedom is apparently left for the numerical parameters at our disposal. Let us 
comment them briefiy. First, the propagator imaginary part S/ has been introduced for 
computational convenience. Hence it can be kept constant when changing the gauge group. 
Second, it can be checked that, by dimensional analysis, our results can all be expressed in 
terms of the ratios T/Tc, T^/ y/a and mo/^/a. 

According to glueball gas models with a Hagedorn spectrum describing the high-lying 



glueball states, the critical temperature is given by |l6l. |43| 

Te [3 




^ - — = 0.69. (53) 

This value is due to the Hagedorn spectrum, not defined above a certain temperature. 
This temperature is here interpreted as the deconfinement one. It is worth saying that it leads 



to an EoS in very good agreement with lattice results 16|, |43|, below Tc. In this picture. 



the ratio T^j \fo is gauge-group independent: This is only valid in a first approximation 
since, for example, there are lattice evidences showing that T(.j ^Ja is only constant up to 
corrections |45|]. Nevertheless, such deviation are beyond the scope of this exploratory 
work. Note that according to ^|, the critical temperature is found to be pretty close to 
300 MeV up to fiuctuation of about 10% for the gauge groups SU(A^), Sp(2), and E7. So, 
we fixed Tc = 300 MeV in our calculations. This value is in good agreement with ( 153|) for 
the value a = 0.176 GeV^ chosen for T = calculations (see Sec. IV BD . 

Concerning the ratio m^j^^ it is worth mentioning the work 47|, in which it is shown 
that the nonperturbative gluon propagator at zero temperature (thus mo in particular) shows 
no significant quantitative differences when expressed in units of the string tension for the 
groups SU(A^) and G2. It is thus tempting to say that the ratio rriQ/ y/a may be gauge-group 
independent also: This is assumed in the rest of this paper. The value mo/-\/o" = 1.67 
obtained from the zero-temperature glueball spectrum is retained. 

Let us finally mention that, at zero temperature, the string breaking scale is found to be 
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two times the gluelump mass for SU(A^) and G2 [18|. This means that the extension of the 
above zero-temperature calculations to any gauge group is straightforward in our approach: 
The T-matrix masses given in Table [TTl can be considered as valid for any gauge group once 
divided by i/a. 



VI. EXISTENCE OF GLUEBALLS ABOVE 



A. Singlet states 

Now that the parameters of the model have been all fixed, T-matrix calculations above 
Tc can be performed. Technical details will not be given here since the method is identical 
to the one used in jsj, which involves the Haftel-Tabakin algorithm to solve the T-matrix 



Lippmann-Schwinger equation [15|. The bound and resonant states appear as poles in the 
on-shell T-matrix, or more precisely as zeros of det J-". The corresponding masses are given in 
Table UTTl in the different considered color-singlet channels. Since in this case, kq-, = —1 for 
all gauge groups and since the gluon mass is independent of A^, the masses of the color-singlet 
are the same for SU(A^) for all A^, and for 6*2 • 

Only a few papers devoted to the existence of glueballs on the lattice are currently 
known 48|], and the interpretation of their results depends mostly on the way the glueball 
correlators are fitted: Either a single narrow pole, or a Breit-Wigner shape. Let us focus 
on the narrow pole fit, which identifies bound states in a way similar to ours. The main 



observation to be made from [48| is that the glueball masses decrease above Tc with increasing 
temperature, with a mass near Tc that is similar to the zero temperature one. This nontrivial 
behavior is well-checked within our approach. Two competing effects are responsible for the 
temperature evolution of the spectrum: reduction of the binding energy and downward shift 
of the threshold energy. Overall, the singlet scalar bound state experiences a mild shift to 
lower energies and dissociates at Tjis ~ 1.3 T^. This is the value from which det J-" does 
not vanish anymore. Nevertheless, considerable strength remains at threshold up to about 
1.5 Tc. This is in qualitative agreement with the spectral function analysis of Euclidean 



correlators by the CLQCD Collaboration [4^ 

The evolution of the imaginary part of the on-shell T-matrix in the singlet scalar channel 
versus the temperature is displayed in Fig. [TJ This gives an overall picture of the glueball 
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TABLE III. Masses (in units of y/a) of lowest-lying glueballs above (T^ = 0.3 GeV). A line mark 
the temperature at which a bound state is not detected anymore. 



Channel 


S 


inglet 


Adjoint-^ 


(2,0)^ 


Group 


All 


SU(Af > 3) 


G2 


T/T, 


2 nig 


0++ 


0-+ 


2++ 


0++ 


0-+ 


2++ 


0++ 


0-+ 


2++ 


1.05 


6.50 


4.52 
6.481 


5.43 


5.43 


6.00 


6.45 


6.31 


6.14 




6.38 


1.10 


5.24 


4.57 


5.21 


5.00 


5.14 






5.21 






1.15 


4.71 


4.43 




4.67 














1.20 


4.43 


4.33 


















1.25 


4.26 


4.24 


















1.30 


4.14 




















1.35 























Radial excitation below the threshold 

progressive dissolution in the medium. The peak in the imaginary part, depicting a bound 
state, becomes broader and broader before melting into the continuum (and thus det J-" does 
not vanish anymore below threshold) as the temperature is increased. Still, for T > Tjis 
and above the threshold energy one finds sizable strength from the bound state relic, the 
T— matrix exhibiting a resonant behavior well beyond the Born approximation. 

Concerning the pseudoscalar channel, singlet bound states are found up to 1.10 Tc. Note 
that states in the pseudoscalar channels, which in our approach correspond to pure P-wave 
scattering, are just mildly bound due to the centrifugal barrier. The tensor states, having 
an S'-wave component, lie between the scalar and pseudoscalar channels, regarding binding 
and dissociation temperatures. 

B. Colored states 

Bound states in the symmetric adjoint channel of SU(A^ > 3) are also observed (see 
Table UTTj) . although they are less bound since K,c-m = The scalar channel disappears 
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FIG. 1. ImT for gg scattering in the 0"''"'' singlet channel for various T with Tc = 0.3 GeV. 



above 1.10 whereas in the pseudoscalar and tensor channels, bound states are lying 
right below the threshold energy at the lowest considered temperature {i.e. 1.05 T^). The 
differences between singlet and adjoint channels have to be attributed to the strength of the 
potential, which is two times smaller in the adjoint channel than in the color singlet. 

The evolution of the T-matrix in the singlet and symmetric adjoint scalar channel versus 
the temperature is displayed in Fig. |2l One clearly sees the disappearance of this bound 
state at 1.15 Tc while the singlet state is still well bound at this temperature. 

There are in general other colored channels than the adjoint one. For SU(A^) gauge groups 
in particular, the only one that could a priori lead to bound states is the (0, 1, 0, . . . , 0, 1, 0)'^ 
channel, which is weakly attractive and exists only for > 3. It has been checked that 
even the scalar state (the most attractive channel) is unbound at = 4. Hence this color 
channel does not admit bound states. In the case of G2, that is the other group considered 
in this study, the (2, 0)'^ channel leads to bound states with the same melting temperatures 
as in the adjoint channel, up to our current precision of 0.05 Tc (see Table UTTll . 
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FIG. 2. T-matrix for gg scattering in the scalar singlet and symmetric adjoint channels for 
SU(A^ > 3). From left to right the temperatures are (1.05;1.10;1.15) T^, with = 0.3 GeV. 



VII. EQUATION OF STATE 

The pure gauge EoS can now be computed without introducing extra parameters: The 
two-body potential and the thermal mass contribution to the gluon mass have been fitted 
on lattice data by using respectively the scaling (146|) and ( l48l) . and they ensure a correct 
agreement with zero temperature results. A crucial point to establish the EoS thanks to f lT6|) 
is to correctly express the summation on the different channels. For 2-gluon interactions, 
the channels are cued by the number and the color number. The summation on the 
channels is formally infinite but in this work, only the ^ and 2++ channels are taken 
into account. This restriction is supported by the following argument. These three channels 
are the most attractive ones and generate the lightest glueballs. The lighter the mass, the 
more important the thermodynamic contribution is in the bound state sector. Thus, the 
bound state thermodynamic contribution coming from other should be negligible in com- 
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parison of these three channels. In order to be coherent, this restriction is implemented 
in the scattering sector. It creates also a restriction in the allowed color channels. Since the 
O"*"^, and 2^^ are symmetric channels, the color channels must be symmetric too in 
order to respect the Pauli's principle. 




TIT, 

FIG. 3. Normalized pressure versus temperature in units of Tc (with Tc = 0.3 GeV), computed for 
the gauge group SU(3) in the free gluon gas case and in the full approach. 

In Fig. [HI the normalized pressure P/Psb computed in the free gluon g with the 

thermal mass ( H9l) is compared with the normalized pressure obtained by our approach in 
the SU(3) case. At low temperature (T < 1.3 Tc), the bound state and the scattering parts 
both give thermodynamic contributions that modify the free gas pressure. For T > l.STc, 
only the scattering part keeps to contribute. As it can be observed in Fig. |3l the main global 
effect of the interaction is to decrease the pressure. If each contribution is analyzed, it is 
seen that the bound state formation increases the pressure because bound states are simply 
added as new species that does not interact with the other particles inside the plasma. 
Concerning the two-gluon scattering part, the sign of the pressure contribution can not 
be analytically predicted at each temperature. Only at the Born approximation, one can 
observe that attractive (repulsive) channels contribute to increase (decrease) the pressure. 
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Indeed, in momentum space representation (see Appendix [C|) . (fTTi) becomes here 



1 foo 2 

^^ = ^4^E(2'^+l)Edi^^'^c,.. / dee'J--mlK,{Pe)vj. (54) 



JP C,g 

where Ki{x) is the modified Bessel function of the first kind. In the attractive (repulsive) 
channels, the sign of the potential is negative (positive). Since Qg is the scattering contribu- 
tion to the grand potential, it can be deduced that attractive (repulsive) channels increase 
(decrease) the pressure. In the present SU(3) case, the only repulsive channel is the (2,2)'^. 
That means that the decreasing of the pressure in our approach compared with the free gas 
pressure is only driven by the (2, 2)'^ channel. 

It is also worth wondering whether some constraints arise or not from the high-temperature 
limit of our framework concerning the behavior of the two-body interactions. In this limit, 
the Born approximation should be relevant. Using fl39|) and i^^, one can write 



Only the scalar channel has been taken into account for the sake of clarity, but the following 
argument can be extended to any spin. According to Hard-Thermal-Loop results, it is 



relevant to assume a Yukawa form for the potential vq at high temperature 



Then, 



^ ^ ^dimo^ r r-^-^'^l .5g. 

where M is the screening mass of the theory, a priori temperature-dependent. Still in HTL 
theory, it is found that, because of the running of the strong coupling constant, 

lim/3^o/3"^9 = lim/3_^o/3M = 0. (57) 

More precisely, the quark and gluon thermal masses are found to behave as as{T)T. Con- 
sequently, at high enough temperatures, it is found that 

- (58) 

i.e. a scattering contribution that has the same behavior with respect to the temperature 
as the free part, ensuring a well-defined large-temperature limit. 

Notice that our fit of the screening mass does not follow the constraints (1571) . but it is 
designed to fit the static potential below 3 T^. A more involved form would be needed to 
reach the HTL predictions at high temperatures, but it is not the scope of the present work. 
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FIG. 4. Normalized pressure versus temperature in units of Tc (with Tc = 0.3 GeV), computed for 
the gauge groups SU(2,3,oo) and G2 (sohd hues). Note that ah the curves are nearly indistinguish- 



able. Our results are compared to the lattice data of 
(dots), and of the minimal G2 model of 



50i] for SU(2) (dots) and 



51| for SU(3,4,6,8) 



52( 1 for G2 (dashed line). Note that all Lattice data have 



been normalized to the lattice Stefan-Boltzmann pressure 



50, 



In Fig. m the normalized pressure P/Psb is presented for different gauge groups: SU(2), 
SU(3), SU(oo) and G2. Severals remarks can be done. First, the free gluon thermodynamic 
contribution is gauge-group invariant once normalized to psB- The gauge-group dependence 
is only present in the bound state and scattering sectors. The number of allowed color 
channels {i.e. the symmetric ones) depends on the gauge group (see ( 1M|) ) and determines 
the allowed maximum number of bound states and the number of scattering channels. The 
bound state thermodynamic contribution comes from two effets: The number and the mass 
of the existing glueballs. Because of the glueball dissociation, this contribution is only 
taken into account up to the temperature of dissociation (see Tables UTTl) . One can observe 
on Fig. m that the produced EoS are not very sensitive to the gauge-group. The most 
important difference between the curves occurs between 1.05 and 1.35 Tc (see Fig. |5]): In 
this range, the gluon-gluon interactions are maximal. When the temperature increases, the 
Born approximation becomes more and more valid and the pressure then scales as dim adj. 
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FIG. 5. Normalized pressure versus temperature in units of Tc (with Tc = 0.3 GeV), computed 
for the gauge groups SU(2), SU(3), G2, and SU(oo). The temperature range is the one where the 
differences between the curves are the most important. 



Thus the normalized pressure tends to be universal. 

In Fig m it is also worth noticing that the EoS computed in our approach favorably 
compares with QCD lattice data for gauge groups SU(3-8) 5l| where such universal curves 
seem t o ap pear (note that lattice data exist also for very high values of T/Tc but only for 
SU(3) 531] )• Concerning G2, no lattice data about EoS are currently available but a new 
effective matrix model describing pure Yang-Mills thermodynamics has been proposed in 



52l |. These last results are compared to ours in Fig HI 



VIII. CONCLUSIONS 



The relevance of gluon-gluon interactions beyond the critical temperature in the pure 
gauge SU(3) plasma has been addressed in a non-perturbative T-matrix many-body frame- 
work with the input of Casimir-scaled potentials from thermal lattice QCD and a model 
of quasigluon mass independent of the gauge group. Scalar glueball bound states in the 
singlet channel survive up to temperatures of about 1.3-1.5 together with sizable thresh- 
old effects due to strong correlations beyond the two-particle threshold. With only one free 
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parameter, the gluon mass at T = 0, the EoS of the gluon-glueball gas is reproduced in good 
agreement with quenched lattice SU(A^) simulations (the other parameters can be fixed a 
priori by resorting to either lattice results or theoretical arguments). Predictions for the G2 
EoS are also given, the main feature being that it should be very close to the SU(N) one 
once normalized to the Stefan-Boltzmann pressure. 

The present T-matrix formalism can in principle be systematically improved, for example 
by including three-gluon scattering. The number of channels under consideration is quite 
larger compared with the two-gluon case, with more elaborated symmetries. Moreover, 
finding the T-matrix would then become a three-body problem, whose resolution through 
e.g. Faddeev equations can be addressed in future works. 

Another natural extension to this paper is to study the light meson spectrum at finite 
temperature and the QCD EoS by including quarks within the model. Computations with 
baryonic potential can be also considered. 

Finally, the T-matrix formalism can also be applied to calculate bulk thermodynami- 
cal properties of the system such as the sheer viscosity, which can be easily computed in 
relaxation-time approximation within a quasi-particle picture. Such a work is in progress. 
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Appendix A: Some SU(A^) relations 

The dimensions of the color channels appearing in ([31]) are given in Table IIV[ together 
with the color factors (l33ll . Note that a general method for computing the quadratic Casimir 



of SU(A^) can be found in 54 1. 



Similar results can be written in the case where only the fundamental and/or conjugate 
representations are taken into account. The terms appearing in the tensor product of the 
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TABLE IV. Symmetry, dimension (dim C) and color factor (kc) defined in (j33p . of the color 
channels (C) appearing in the tensor product of the SU(A^) adjoint representation by itself. This 
table actually displays the two-gluons color channels, denoted by gg. The SU(3) case is also 
indicated. 

Ciorgg • (1, 0, . . . , 0, 1) (2, 0, . . . , 0, 2) (2, 0, . . . , 1, 0) (0, 1, 0, . . . , 0, 1, 0) 

(0,1,..., 2) 

SU(3) . (1,1) (2,2) (0,3), (3,0) 



Symmetry 
dim C 



S 
1 

-1 



S, A 
_ 1 

2 



S A 

N^{N+3){N~l) {N^-A){N^-1) 
4 4 



1_ 
N 







' N 



fundamental and conjugate representations by themselves are given in Table |Vl as well as 
those appearing in the tensor product of the fundamental representation by the conjugate 
one. 



TABLE V. Symmetry, dimension (dim C) and color factor (kq) defined in ([8|), of the color channels 
(C) appearing in the tensor product of the SU(A^) fundamental representation by itself (left), of the 
conjugate representation by itself (middle), and of the fundamental representation by the conjugate 
one (right). This table actually displays the quark-quark, antiquark-antiquark and quark-antiquark 
color channels, denoted by qq, qq, and qq respectively. 



C for qq, qq, qq 


(2,0,..., 


0) (0,1,0,... 


0) 


(0,...,0,2) (0, 


...,0,1,0) 


• 


(1,0,... ,0,1) 


Symmetry 


S 


A 




S 


A 






dim C 


N{N+1) 
2 


N{N~l) 
2 




N{N+1) 
2 


N{N-1) 
2 


1 


iV2-l 




N-l 


N+1 




N-l 


N+l 


V2-l 


1 


2N'^ 


2N'^ 




2N'^ 


2N'^ 


2N'^ 


2N^ 



Finally, useful results concerning the tensor product of the fundamental (conjugate) rep- 
resentation by the adjoint representation are given in Table IVII 



Appendix B: Lattice potential 



The lattice data that are used as a starting point to build our interaction potential are 
those of (22|, i.e. the static free energy between a quark-antiquark pair bound in a color 
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TABLE VI. Dimension (dim C) and color factor (kc) defined in ^ of the color channels (C) 
appearing in the tensor product of the SU(A^) fundamental representation by the adjoint one 
(left), and of the conjugate representation by the adjoint one (right). This table actually displays 
the quark-gluon and antiquark-gluon color channels, denoted by qg and qg respectively. 



C for qg, qg 


(1,0,... 


,0) (2,0,... ,0,1) (0,1,. ..,0,1) 


(0,. 


..,0,1) 


(1,0,... ,0,2) (1,0,...1,0) 


dim C 


N 


{N+2)N{N-l) {N+l)N(N-2) 




N 


(Af+2)Af(Af-l) (Af+l)Af(Af-2) 


2 2 




2 2 




1 

2 


1 1 

2N 2N 




1 

2 


1 1 

2N 2N 



singlet for Niat = 3. For numerical convenience, it is preferable to deal with a fitted form 



of these, rather than with interpolations of the available points. To fit the data of 



22], the 



analytic form proposed by Satz in 55[ is used: 

r(i/4) VmtT 



Fi(r,T) 



a 



/i(T) 



(^(T)V2) 



4a 
3 r 



-M(T)r 



+ /i(T)r . (Bl) 



23/2r(3/4) 23/4r(3/4)' 

The way of obtaining this formula is the following. First, it is known that the static quark- 
antiquark energy at zero temperature is accurately fitted by a so-called funnel shape 

4a 



Fi(r, 0) = ar — 



3r 



f/i(r,0), 



(B2) 



see e.g. [37[. When T > 0, one can imagine that this potential is progressively screened by 
thermal fluctuations. An effective theory for studying the screening of a given potential is 
the Debye-Hiickel theory, in which the thermal fluctuations are all contained in a screening 
function /i(T), that modifies the zero-temperature potential and eventually leads to the form 

m- 

The explicit form of /i(T) is unknown a priori and has to be fitted on the lattice data. 
As it can be seen in Fig. |6l the form 



/i(T) 



T 



0.537— + 0.644 + 0.112 In — - 0.967 



T 



(B3) 



with 



a 



0.141, (B4) 

provides an accurate fit of the lattice data in the range 1-3 Tc. A more complete fit should be 
such that /i(0) = 0, but our model is not intended to be able to "cross" the phase transition 
in Tc. The simple form (]B3|) is already satisfactory. The corresponding internal energy 
Ui = Fi- TdrFi is plotted in Fig. El 
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FIG. 6. Static free energy Fi{r,T) of a quark-antiquark pair bound in a color singlet, computed 
in SUf3) quenched lattice QCD and plotted for different temperatures (symbols). Data are taken 



from 



22| and expressed in units of with r the quark-antiquark separation. The fitted form 



()Bip - (|B4p is compared to the lattice data (solid lines). 




Appendix C: Dashen, Ma and Bernstein's formalism in momentum space 



To compute ( |T6l) . it is necessary to use a given representation. In order to use the 
calculation of the T-matrix proposed in [sl, the scattering part of ( ITTj) must be computed 
in the momentum space representation (the two first term are simply free gas contributions 
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and can be easily computed). Let us focus on 

dimC 



C JP 

xTrcjp 



(SReT)' - 2n({5Rer)i5lmry - {5lmr){5Rer)' 



(CI) 



Using the following definitions concerning the trace of an operator A in momentum space 

1 



TiA 



{2nY 



dq {q\A\q) 



and the partial wave expansion 

{q\A\q') = A{q, q\ q.q') = ^ ^^(2/ + l)A{q, q')Pi{q.q' 

I 

where Pi{x) is the Legendre polynomial of order /, ( ICip reads 



(C2) 



(C3) 



647r 



+ 



k^J2^2J+l)Y,dimC L rfe ^ - m2J^i(/3e) ReTjp (e; g„ g,) 

jp c \ ^ 

K2{l3e) [Re7}p(e; g,, g,) (Im7}p(e; g,, g,))'] 

[(Rerjp(e;g„g,))'lmrjp(e;g„g,)] ) . (C4) 



^ , dee' \ m„ 

167r2 I 4 ^ 



^ / ae e - — m„ 
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